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ABSTRACT 

This study focuses on some of the most metal-poor damped Lya absorbers known 
in the spectra of high redshift QSOs, using new and archival observations obtained with 
UV-sensitive echelle spectrographs on the Keck and VLT telescopes. The weakness 
and simple velocity structure of the absorption lines in these systems allow us to 
measure the abundances of several elements, and in particular those of C, N, and 
O, a group that is difficult to study in DLAs of more typical metallicities. We find 
that when the oxygen abundance is less than ~ 1/100 of solar, the C/O ratio in 
high redshift DLAs and sub-DLAs matches that of halo stars of similar metallicity 
and shows higher values than expected from galactic chemical evolution models based 
on conventional stellar yields. Furthermore, there are indications that at these low 
metallicities the N/O ratio may also be above simple expectations and may exhibit a 
minimum value, as proposed by Centurion and her collaborators in 2003. Both results 
can be interpreted as evidence for enhanced production of C and N by massive stars in 
the first few episodes of star formation, in our Galaxy and in the distant proto-galaxies 
seen as QSO absorbers. The higher stellar yields implied may have an origin in stellar 
rotation which promotes mixing in the stars' interiors, as considered in some recent 
model calculations. We briefly discuss the relevance of these results to current ideas 
on the origin of metals in the intergalactic medium and the universality of the stellar 
initial mass function. 
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1 INTRODUCTION 

The neutral hydrogen clouds which give rise to damped Ly- 
man alpha systems (DLAs; Wolfe, Gawiser, & Prochaska 
2005) in the spectra of high redshift QSOs provide impor- 
tant information on the early chemical evolution of galaxies, 
complementing that obtained from more traditional studies 
of metal-poor stars and Hn regions in the Milky Way and 
nearby galaxies (Pettini 2004; Henry & Prochaska 2007). 
The relatively simple physics of interstellar line formation 
(Stromgren 1948) makes the determination of ion column 

* Based in part on data obtained at the W.M. Keck Observa- 
tory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of Califor- 
nia, and NASA, and was made possible by the generous finan- 
cial support of the W.M. Keck Foundation. Based in part on 
UVES observations made with the European Southern Observa- 
tory VLT/Kueyen telescope at Paranal, Chile, obtained in pro- 
gramme 078.A-0185(A) and from the public data archive. 



densities and corresponding element abundances straight- 
forward, and the rich variety of resonant transitions at rest- 
frame ultraviolet wavelengths allows a detailed picture of 
the chemical composition of many DLAs to be assembled 
(e.g. Prochaska, Howk, & Wolfe 2003; Rix et al. 2007). The 
elements carbon, nitrogen and oxygen have, however, been 
comparatively little studied so far in DLAs, despite the key 
role they play in stellar nucleosynthesis. C and O are abun- 
dant elements with intrinsically strong atomic transitions; in 
combination, these two factors result in resonance absorp- 
tion lines which are strongly saturated — and thus unusable 
for abundance determinations — even at the low metallici- 
ties of most DLAs. The problem is particularly acute for 
C. For O some weaker transitions are available at far-UV 
wavelengths; these, however, can only be accessed from the 
ground at redshifts z ^ 2.5. The absorption lines of N, on 
the other hand, tend to be weak and are often blended with 
intergalactic Lya forest lines; nevertheless, the abundance of 
N in DLAs has been the subject of several studies since at- 
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Table 1. Journal of Observations 



QSO 


V or g a 






Telescope/ 


Wavelength 


Resolution 


Integration Time 


S/N b 


7V(H i) 




(mag) 






Instrument 


Range (A) 


(km s _1 ) 


w 




(cm~ 2 ) 


Q0913+072 


17.1 


2.785 


2.61843 


VLT/UVES 


3310-9280 c 


6.7 


77 500 


60 


2.2 X 10 20 


SDSS J1016+4040 


19.5 


2.991 


2.81633 


Kecki/HIRESb 


3190-5970 c 


7.5 


15 000 


12 


0.8 x 10 20 


SDSS J1558+4053 


18.7 


2.635 


2.55332 


Kecki/HIRESb 


3190-5970 c 


7.5 


25 200 


18 


2.0 x 10 20 


Q2206-199 


17.3 


2.559 


2.07624 


VLT/UVES & 


3100-10 000 c 


6.5 


68 500 


100 


2.7 x 10 20 










Kecki/HIRESr 













a Magnitudes are g for SDSS objects and V for other QSOs. 
b Indicative signal-to-noise ratio at 5000 A. 
c With some wavelength gaps 



tention was first drawn to it by Pettini, Lipman & Hunstead 
(1995). 

The relative abundances of C, N, and O at low metal- 
licities hold clues to the nature of the stars responsible for 
their production in the earliest stages of nucleosynthesis in 
galaxies, and their measurement in DLAs has the potential 
of clarifying some unresolved questions raised by recent stel- 
lar work, as discussed in detail later. The marked increase 
in the number of known DLAs afforded by the Sloan Digital 
Sky Survey (Prochaska, Herbert-Fort, & Wolfe 2005) makes 
it possible to search for damped systems with the character- 
istics which facilitate the measurement of the abundances of 
C, N, and O, primarily low metallicity and simple velocity 
structure. In this paper we report new high resolution spec- 
troscopic observations of four such absorption systemsQ all 
with metallicities [O/H] < —20 When combined with pub- 
lished measurements in two similar systems, this data set — 
although modest — sheds new light on the nucleosynthesis of 
C, N, and O in the low metallicity regime. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target Selection 

The four QSOs observed were selected based on the weak- 
ness of metal lines in DLAs already known to be present in 
their spectra, either from the trawl of the SDSS QSO sample 
by Prochaska et al. (2005) or from earlier published stud- 
ies. It has now become apparent that metallicity and strong 
metal line widths are correlated in DLAs (e.g. Ledoux et al. 
2006a; Murphy et al. 2007a; Prochaska et al. 2007). Thus, an 
efficient way to preselect for further study the most metal- 
poor DLAs is to focus on DLAs in which the associated 
strongest metal lines, which are nearly always saturated, 
have very low values of equivalent width — so low as to re- 
main undetected at the moderate resolution and signal-to- 
noise ratio (S/N) of most SDSS spectra. In addition, the 
simple kinematic structure of the absorbing gas implied by 
the low equivalent widths simplifies the abundance analysis 
compared with the more usual occurrence of multiple ab- 
sorption components spread over 100-200 km s" 1 . 

1 Three DLAs and one sub-DLA, according to the conventional 
distinction between the two at neutral hydrogen column densities 
JV(Hl)= 2 x 10 20 cm- 2 . 

2 As usual, [O/H] = log(0/H) DLA - log(O/H) . 



Using this method, we have discovered two of the 
most metal-poor DLAs, with oxygen abundance [O/H] ~ 
-2.4, in the spectra of the QSOs SDSS J1016+4040 and 
SDSS J1558+4053; high resolution observations of these sys- 
tems are reported here for the first time. We have also re- 
observed two QSOs with known low-metallicity, narrow-line 
DLAs: Q0913+072 and Q2206-199; the new data add sig- 
nificantly to those already available in the literature thanks 
to the very high S/N ratio achieved. 

2.2 Observations 

The observations were performed with the HIRES (Vogt et 
al. 1994) and UVES (Dekker et al. 2000) echelle spectro- 
graphs on, respectively, the Keck I and VLT2 telescopes. Rel- 
evant details are collected in Table Q] For the SDSS QSOs 
we used the refurbished HIRES with its UV-sensitive cross- 
disperser, covering the wavelength range 3190-5970 A. The 
UVES observations of Q0913+072 and Q2206-199 used a 
variety of dichroics and wavelength settings to cover essen- 
tially all of the optical wavelength range from ~ 3200 — 
3300 A to ~ 9000 - 10 000 A. The resolution is FWHM ~ 
6 — 7 km s _1 for spectra recorded with either spectrograph. 
The S/N ratio varies considerably along the spectrum, par- 
ticularly for the UVES data where a given portion of the 
spectrum may have been covered with a number of different 
instrument settings; in the penultimate column of Table [1] 
we list an indicative value of S/N at 5000 A. 

2.3 Data reduction 

The reduction of the two dimensional images into co-added 
one-dimensional spectra followed standard procedures. The 
HIRES spectra were reduced with the makee pipeline soft- 
ware kindly maintained by T.A. Barlow, while for the UVES 
spectra we used the pipeline reduction provided by the Eu- 
ropean Southern Observatory but with modifications which 
improve the accuracy of the wavelength calibration and flux 
extraction (see Murphy et al. 2007b). These pipelines take 
care of the initial stages of bias subtraction, flat-fielding, 
wavelength calibration, extraction of individual echelle or- 
ders, and flux calibration. The individual spectra of each 
QSO were averaged and the echelle orders merged with 
purpose-built software^ designed to maximise the signal-to- 

3 Available from http: / / astronomy.swin.edu.au/ ~mmurphy/UVES_popler 
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Figure 1. Observed profiles (black histograms) and fitted Voigt profiles (continuous red lines) of selected absorption lines in the 
•Jabs = 2.6f843 DLA in Q0913+072. Dashed red lines indicate nearby absorption not attributed to the DLA. The y-axes of the plots 
show residual intensity; note the expanded y-scale of the bottom plot. The full list of absorption lines analysed in this DLA is given in 
Table [3] 
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noise ratio of the final spectrum which was binned onto a 
linear, vacuum heliocentric wavelength scale with approxi- 
mately three bins per resolution element. 

2.4 Abundance Measurements 

Absorption lines from each DLA were identified and their 
profiles fitted with Voigt profiles generated by the VPFIT 
(version 8.02) software package, as described by Pettini et al. 
(2002a) and Rix et al. (2007). Briefly, VPFIT uses x" 2 min- 
imisation to deduce the values of redshift z, column density 
N (cm -2 ), and Doppler parameter b (km s _1 ) that best re- 
produce the observed absorption line profiles. VPFIT takes 
into account the instrumental broadening function in its 
X 2 minimisation and error evaluation^ We used the com- 
pilation of laboratory wavelengths and /-values by Morton 
(2003) with recent updates by Jenkins & Tripp (2006). 

The DLAs observed are among those with the simplest 
kinematic structure, generally consisting of only two or three 
absorption components separated by ~ 10 — 20 km s — 1 ; 
this fact, coupled with the high resolution and generally 
high S/N of the spectra, made the profile fitting procedure 
straightforward and its results unambiguous. Another conse- 
quence of the unusually simple line profiles is that the same 
set of 'cloud' parameters, b and z, was found to reproduce 
well all the absorption lines considered in each DLA; nat- 
urally, lines arising from the same ground state of a given 
ion were fitted together to give the best-fitting value of the 
column density N(X 1 ), where N(X l ) is the ith ion stage 
of element X. In each DLA, the column density of neutral 
hydrogen was determined by fitting the damping wings of 
the Lya line which are its most sensitive measure; we also 
checked that the resulting value of 7V(H i) provides a good 
fit to the higher order Lyman lines covered by our spectra. 

Element abundances were deduced from the observed 
ratios iV(X n)/N(H i), where N(Xn) denotes the column 
density of the ion stage which contains the dominant frac- 
tion of element X in H I regions of the Milky Way interstellar 
medium (e.g. Morton et al. 1973). For the elements consid- 
ered here, these are either the neutrals (Ol, Nl) or the first 
ions (C II, Alii, Sill and Fell). However, we have also mea- 
sured the column densities of higher ions stages, when their 
absorption lines are present in our spectra, in order to as- 
sess the corrections, if any, that need to be applied to our 
abundance determinations to account for the possible con- 
tributions by ionised gas to the absorption lines observed. 
When element abundances are referred to the solar scale, we 
have used the compilation by Asplund, Grevesse, & Sauval 
(2005). 



3 INDIVIDUAL ABSORPTION SYSTEMS 

3.1 Q0913+072; DLA at z abs = 2.61843 

The metal-poor DLA in front of this bright QSO has been 
known since the survey by Pettini et al. (1997). It was sub- 
sequently studied at higher spectral resolution by Ledoux et 
al. (1998) and Erni et al. (2006) who considered it to be 'the 

4 VPFIT is available from http://www.ast.cam.ac.uk/"' 
rfc/vpfit.html 



Table 2. Absorption Components of Low Ion Transitions 
in Q0913+072 



Component 


2 abs 


b 


Fraction a 


Number 




(km s" 1 ) 




1 


2.61828 


3.7 


0.28 


2 


2.61843 


5.4 


0.72 



a Fraction of the total column density of Si II. 



Table 3. Ion Column Densities in Q0913+072, z ahs = 2.61843 
DLA 



Ion 




Transitions used 


logJV(X) 


Hi 


1215 




20.34 ±0.04 


Cn 


1036, 


1334 


13.98 ±0.05 


N I 


1199 




12.29 ±0.12 


Oi 


1039, 


1302 


14.63 ±0.01 


Alii 


1670 




11.78 ±0.03 


Aim 


1854 




< 11.03 


Sin 


1190, 


1193, 1260, 1304, 1526 


13.30 ±0.01 


Fen 


1608, 2344, 2374, 2382 


12.99 ±0.01 



most metal- deficient DLA known' and drew attention to its 
role as a probe of early nucleosynthesis. As part of a pro- 
gramme to measure the primordial abundance of deuterium, 
we re-observed this QSO with UVES in 2007. Including a 
few spectra downloaded from the UVES data archive, the 
total exposure time devoted to Q0913+072 is 77 500 s; the 
corresponding signal-to-noise ratio near 5000 A is S/N~ 60 
(see Table 1). 

A selection of the absorption lines in this DLA is shown 
in Figure [T] The metal lines consist of two closely spaced 
components; profile decomposition with VPFIT returns b- 
values of 3.7 and 5.4 km s" 1 and a velocity separation of 
12.4km s _1 (see Tabled- Table [3] lists the absorption lines 
analysed and the corresponding values of column density; 
the errors on N(X) are the formal values returned by VP- 
FIT. Most of the metal lines in this DLA are weak; in partic- 
ular, from Figure [T] it can be seen that both CnA1334 and 
C II A1036 are unsaturated (these lines are nearly always sat- 
urated in DLAs), as is OiA1039. Although both Cn lines 
are partially blended with broad absorption at longer wave- 
length (modelled by VPFIT with two components), the im- 
pact of these contaminating features on the determination 
of N(C n) is relatively minor: a conservative estimate of the 
error on log V(C n) deduced from several VPFIT trials is 
±0.05. 

We have a definite (4a) detection of N I Al 199.55 in both 
velocity components; the two other, weaker, members of the 
Nl AA1199.55, 1200.22, 1200.71 triplet are blended with Lya 
forest lines. 

Turning to absorption lines from ionised gas, CmA977, 
Si in A1206, and FeniA1122 are blended with Lya forest 
lines, while NnA1084 is undetected (which is unsurprising 
given the weakness of the Nl lines). We do, however, have 
a tentative detection of AlniA1854; given the weakness of 
this feature, we consider the column density of Aim to be 
an upper limit. 
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Figure 2. Observed profiles (black histograms) and fitted Voigt profiles (continuous red lines) of selected absorption lines in the 
•Jabs = 2.8f633 DLA in SDSS J1016+4040. Dashed red lines indicate nearby absorption not attributed to the DLA. The j/-axes of the 
plots show residual intensity. The full list of absorption lines analysed in this DLA is given in Table [5] 
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Table 4. Absorption Components of Low Ion Transitions 
in SDSS J1016+4040 



Component 


^•abs 


b 


Fraction a 


Number 




(km s _1 ) 




1 


2.81604 


5.0 


0.11 


2 


2.81633 


6.9 


0.89 



a Fraction of the total column density of C II. 



Table 5. Ion Column Densities in SDDS J1016+4040, z abs = 
2.81633 DLA 



Ion 


Transitions used 


log iV (X) 


Hi 


1215 


19.90 ±0.11 


Cii 


1334, 1036 


13.66 ±0.04 


Nl 


1199.5 


< 12.76 


Ol 


1302, 988.77 


14.13 ±0.03 


Sin 


1260, 1526 


12.90 ±0.05 



3.2 SDSS J1016+4040; DLA at z abs = 2.81633 

This QSO is one order of magnitude fainter than 
Q0913±072. We recorded its spectrum from 3190 A to 
5970 A on the night of 2007 May 6 with the refurbished 
HIRES spectrograph on the Keck I telescope, using the UV 
cross-disperser and UV sensitive CCD array. With a total 
exposure time of 15 000 s we achieved a modest S/N~ 12 
near 5000 A. As can be seen from Figure [2] what appears as 
a likely damped Lya line in the low resolution SDSS spec- 
trum breaks up into three distinct components at the higher 
resolution of the HIRES data. However, the central one of 
these three absorption lines does exhibit damping wings; we 
deduce A/(Hi) = (8±2) x 10 19 cm -2 . Such systems are some- 
times referred to as sub-DLAs since their column densities 
are lower than the historical definition of a DLA as a system 
with AT (Hi) > 2 x 10 20 cm" 2 (Wolfe et al. 1986). 

All the metal lines in this DLA are weak (see Figure [5}. 
The first ions exhibit two narrow velocity components with 
b = 5.0 and 6.9 km s _1 separated by 23 km s , with the 
latter accounting for ~ 90% of the total column density 
(see Table Q}. The former, at z abs = 2.81604, probably 
arises in ionised gas since it is not detected in the OlA1302 
and A988.77 lines (the latter being the strongest member of 
the O I AA988.77, 988.65, 988.58 triplet shown in Figure [2]). 
Therefore, we have not included this component in the val- 
ues of N(C II) and N (Si n) attributed to the DLA in Tabled 
Nl is below our detection limit; similarly, absorption lines 
from higher ion stages than those dominant in Hi regions 
are either blended (e.g. SimA1206), or undetected. 

3.3 SDSS J1558+4053; DLA at z abs = 2.55332 

We observed this QSO with HIRES on the night of 2006 
June 17, achieving S/N~ 18 near 5000 A with a 25 200 s ex- 
posure. The echelle data confirm that the strong absorption 
feature near 4320 A in the discovery SDSS spectrum is in- 
deed a damped Lya line with N(H i) = 2.0 x 10 20 cm -2 (see 
Figure [3]). 

The neutral gas is concentrated in a single, narrow 



Table 6. Absorption Components of Low Ion Transitions 
in SDSS J1558±4053 



Component 


■^abs 


b 


Fraction a 


Number 




(km s- 1 ) 




1 


2.55264 


9.1 


0.01 


2 


2.55298 


8.8 


0.06 


3 


2.55332 


6.1 


0.80 


4 


2.55334 


25.1 


0.13 



a Fraction of the total column density of Si II. 



Table 7. Ion Column Densities in SDDS J1558±4053, z abs = 
2.55332 DLA 



Ion 


Transitions used 


logA(X) 


Hi 


1215 


20.30 ±0.04 


Cii 


1334, 1036 


14.22 ±0.06 


Nl 


1199.5, 1200.2, 1200.7 


12.66 ±0.07 


Ol 


1302, 988.77, 988.65 


14.54 ±0.04 


Alii 


1670 


11.92 ±0.06 


Sin 


1190, 1193, 1260, 1304, 1526 


13.32 ±0.02 


Fen 


1096, 1144, 1260, 1608, 


13.07 ±0.06 



(b = 6.1km s _1 ) component at ;z abs = 2.55332; other, 
weaker, components are seen at positive and negative ve- 
locities in C II and Si II and probably originate in ionised 
gas (Figure [3] and Table [6j . The entries in Table [7] are for 
the neutral component only. The absorption lines in this 
DLA are stronger than in the previous two cases. Neverthe- 
less, we cover at least one unsaturated transition for every 
species considered except possibly for Cll: both A1334 and 
A1036 are just saturated. Even for C II, however, we can de- 
termine the column density with confidence because: (i) the 
absorption profile (of the central, dominant component) is 
well fitted by the same b- value which fits best all the other 
lines; (ii) larger values of b are incompatible with the ob- 
served width of the line which is resolved in our data. Note 
that while we cannot exclude, in principle, that narrower 
components may be present within the marginally saturated 
C II line profiles, their unrecognised presence would lead to 
an underestimate of [C/H] and C/O, further strengthening 
the main conclusion of our analysis of the behaviour of the 
C/O ratio at low metallicities, as discussed in Section [5] 

As for the case of the DLA in SDSS J1016±4040, we are 
unable to deduce reliable column densities of higher ionisa- 
tion stages than those listed in Table [7] because the relevant 
absorption lines are either blended, undetected, or outside 
the wavelength range of our HIRES spectrum. 

3.4 Q2206 199; DLA at z abs = 2.07624 

The metal-poor DLA in this QSO has been known since the 
early echelle observations by Pettini et al. (1990) and by 
Prochaska & Wolfe (1997). Here we combine spectra from 
the UVES and HIRES data archives to achieve S/N~ 100 
at 5000 A and extensive spectral coverage from ~ 3100 A to 
1 fim. 

Even at this high resolution and signal-to-noise ratio 
this DLA appears to consist of a single, narrow component: 
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Figure 3. Observed profiles (black histograms) and fitted Voigt profiles (continuous red lines) of selected absorption lines in the 
•Jabs = 2.55332 DLA in SDSS J1558+4053. Dashed red lines indicate nearby absorption not attributed to the DLA. The i/-axes of the 
plots show residual intensity. The full list of absorption lines analysed in this DLA is given in Table [7] 
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Figure 4. Observed profiles (black histograms) and fitted Voigt profiles (continuous red lines) of selected absorption lines in the 
Zshs = 2.07624 DLA in Q2206— f99. Dashed red lines indicate nearby absorption not attributed to the DLA. The y-axes of the plots 
show residual intensity. The full list of absorption lines analysed in this DLA is given in Table [9] 
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Table 8. Absorption Components of Low Ion Transitions 
in Q2206-199 



Component 




b Fraction a 


Number 




(km s- 1 ) 


1 


2.07624 


6.5 1.00 



a Fraction of the total column density of Si II. 



Table 9. Ion Column Densities in Q2206-199, z abs = 2.07624 
DLA 



Ion 


Transitions used 


log Af (X) 


Hi 


1215 


20.43 ±0.04 


Cn 


1334, 1036 


14.41 ±0.03 


N I 


1199.5, 1200.2, 1200.7 


12.79 ±0.05 


Nn 


1083 


12.46 ±0.10 


Oi 


1302, 1039 


15.05 ±0.03 


Alii 


1670 


12.18 ±0.01 


Aim 


1854, 1862 


11.59 ±0.04 


Sin 


1193, 1260, 1304, 1808 


13.65 ±0.01 


Sim 


1206 


< 12.86 


Fen 


1121, 1260, 1608, 2344, 2374, 2382, 2600 


13.33 ±0.01 



the profiles of all the absorption lines from species which 
are dominant in H I regions are well reproduced by a model 
'cloud' with b — 6.5 km s _1 (Figure [4] and Tables [8] and 
[9|. The availability of a number of transitions with differ- 
ing /-values for most ions ensures that the column densities 
are well determined. The Cn lines are just saturated but, 
for the reasons outlined above, this is unlikely to affect our 
conclusions concerning the relative abundances of C and O. 
Thanks to its high S/N ratio and wide wavelength coverage, 
our combined UVES±HIRES spectrum does cover a number 
of transitions from ions that only occur in Hn regions (Nil, 
Al m, and Si ill; see Table|9} ; their column densities are lower 
than those of the corresponding Hi region species. Finally, 
it can be seen from Figure|4]that we have a clear detection of 
all three lines of the N I AA1199.55, 1200.22, 1200.71 triplet. 



4 ELEMENT ABUNDANCES 

Our abundance measurements in the four DLAs are collected 
in Table 1101 We have searched the literature for other ex- 
amples of DLAs (or sub-DLAs) where the C II and O I ab- 
sorption lines are unsaturated, so that the corresponding 
column densities can be determined, and found only two 
other cases, also included in Table 1101 from the work by 
Dessauges-Zavadsky et al. (2003) and Peroux et al. (2007). 
Both of these cases are of somewhat lower hydrogen column 
density and at higher redshifts than the four reported here 
but, given the paucity of C/O measurements up to now, we 
shall consider them in the present analysis. 

As can be seen from Table 1101 all four DLAs in the 
present work have metallicities of less than 1/100 of solar. 
The DLA in front of Q0913±072 remains the most metal- 
poor in [Fe/H], but its oxygen abundance of ~ 1/250 solar 
is similar to those of the two newly discovered DLAs from 
the SDSS. With [O/H] ~ -2.4 and [C/H] 2.5 to -2.7, 



these DLAs have attained a level of metal enrichment com- 
parable to that typical of intergalactic Lya forest clouds at 
the same redshifts (e.g. Simcoe et al. 2004). Although the 
source of these intergalactic metals is still somewhat contro- 
versial (see Ryan- Weber et al. 2006 and references therein), 
it seems likely that such metal-poor DLAs have experienced 
very little, if any, star formation in situ. 

It is also interesting to note that no DLA has yet been 
found with metallicity less than 1/1000 solar, despite the 
fact that such metal-poor stars do exist in the halo of our 
Galaxy. In drawing such comparisons, however, it must be 
borne in mind that such extremely metal-poor stars consti- 
tute a tiny fraction of the stellar populations of the Milky 
Way. It may simply be the case that with only a few hundred 
metallicity measurements in DLAs so far, we simply do not 
have the statistics necessary to pick out such outliers in the 
metallicity distribution. 

Finally, we draw attention to the fact that in the four 
systems in which this ratio could be measured, [O/Fe] ~ 
+0.25 to +0.5 (see Table 1 10 p consistent with the enhance- 
ment of the alpha-capture elements relative to iron seen in 
halo stars of similar metallicity. At these low DLA metal- 
licities, the depletion of Fe onto dust grains is expected to 
be minimal (Vladilo 2004; Akerman et al. 2005) so that the 
observed N(0 l)/N(Fe n) ratio gives a good indication of 
inherent departures from the solar relative abundances of 
these two elements. As we show in section 14.11 differential 
ionisation corrections between O I and Fe II are also unlikely 
to account for this effect. 



4.1 Ionisation Corrections 

As explained in section [2~4l the relative element abundances 
listed in Table [101 were derived by dividing the column den- 
sities A^(Xn) by A r (Hl), where W(Xn) is the dominant ion- 
isation stage of element X in Hi regions. Strictly speak- 
ing, this commonly adopted assumption only holds true 
for O out of the elements considered here. It has been 
known for a long time that charge-exchange reactions keep 
N(Ol)/N(Rl) very close to O/H for a wide range of phys- 
ical conditions (e.g. Field & Steigman 1971). For elements 
which in H I regions are mostly singly ionised — that is, C, Al, 
Si and Fe in our case — the possibility exists that some Hn 
gas — at the same velocity as the H I and thus indistinguish- 
able kinematically — may contribute to column densities of 
their first ions; if true, this would lead to an overestimate 
of the corresponding element abundances^] Conversely, N 
may be overionised relative to H, so that even in H I regions 
N(Nl)/N(Rl) < (N/O) leading to an tmderestimate of the 
N abundance. 

Of course, one of the major motivations for using DLAs 
in abundance studies, is the fact that when the column den- 
sity of neutral gas exceeds A r (Hl) = 1 x 10 20 cm -2 such ion- 
isation corrections are expected to be small and comparable 
to other sources of error, as concluded by a number of de- 
tailed analyses (e.g. Viegas 1995; Vladilo et al. 2001). Given 

5 It is worthwhile pointing out here that the neutral atoms of 
these species, which we have neglected, make negligible contribu- 
tions to the abundance measurements, because even in H I regions 
7V(C i) -C N(C II), and similar inequalities apply to Al, Si and Fe. 
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Table 10. Element Abundances 



QSO 




logA^H i) 


[0/H] a 


log (C/0) b 


log (N/0) c 


[Al/H] a 


[Si/H] a 


[Fe/H] a 


Q0913+072 


2.61843 


20.34 


-2.37 


-0.65 


-2.34 


-2.99 


-2.55 


-2.80 


SDSS J1016+4040 


2.81633 


19.90 


-2.43 


-0.47 


sC -1.37 




-2.51 




SDSS J1558+4053 


2.55332 


20.30 


-2.42 


-0.32 


-1.88 


-2.82 


-2.49 


-2.68 


Q2206-199 


2.07624 


20.43 


-2.04 


-0.64 


-2.26 


-2.69 


-2.29 


-2.55 


Other Objects: 


















J0137-4224 d 


3.665 


19.11 


-2.39 


-0.24 






-2.22 




J2155+1358 c 


4.21244 


19.61 


-1.77 


-0.55 




-2.12 


-1.87 


-2.13 



a [X/H] = log (X/H) DLA - log (X/H) Q . Solar abundances from Asplund et al. (2005). 

b These values are not referred to the solar abundance. For reference, log (C/0)q = —0.27 (Asplund et al. 2005). 
c These values are not referred to the solar abundance. For reference, log (N/0)g = —0.88 (Asplund et al. 2005). 
d Peroux et al. (2007). 
°Dessauges-Zavadsky et al. (2003). 



the importance of this point for the measurement of the C/O 
ratio, however, we verified that this is indeed the case for the 
two DLAs in our small sample where column densities of ions 
from H II regions could be measured in our spectra. 

Specifically, we run a series of photoionization models 
with the software package CLOUDY (Ferland et al. 1998; Fer- 
land 2000 Jf], approximating the DLAs as slabs of constant 
density gas irradiated by the metagalactic ionising back- 
ground (Haardt & Madau 2001) and the cosmic microwave 
background. We assumed a metallicity of 1/100 solar and 
no differential depletions due to dust grains (Akerman et 
al. 2005) , but the results of the computations are not sensi- 
tive to these assumptions. We varied the gas volume density 
from log[n(H)/cm" 3 ] = —3 to +3 and stopped the calcu- 
lation when the measured column density of neutral gas in 
each DLA was reached. From the CLOUDY output we then 
calculated for each element observed the ionization correc- 
tion 



IC(X) = log 



JV(X) 



JV(H) 



log 



AT(Xn) 



iV(Hi) 



(1) 



computed 

which is negative when the abundance of an element has 
been overestimated by neglecting the ionized gas. 

Figure [S] illustrates the results of this exercise for the 
DLA in Q2206— 199 (the conclusions are similar for the 
other DLAs considered here). As expected, Ol tracks Hi 
very closely at all densities. At very low densities N is ove- 
rionised relative to H (that is some N is present as N II even 
in Hi regions), while the column densities of the other ele- 
ments considered have some contribution from ionised gas so 
that, for example, N(C li)/N(H i) > C/H. However, for the 
elements which are the focus of the present work, C, N (and 
O), the ionisation corrections amount to less than 0.1 dex 
provided the gas density is greater than n(H) ~ 0.01 cm" 3 . 
When n(H) ^0.1 cm" 3 , the ionisation corrections for all the 
elements in Table 1101 are 0.05 dex. 

We can obtain an indication of the density which ap- 
plies to the gas in the DLAs from consideration of the ratios 
of different ions stages, since they also depend on the den- 
sity, once the radiation field has been fixed (see Figure [6]). 

6 http://www.nublado.org/ 



The DLA which provides the most extensive information 
in this respect is the one in line to Q2206— 199. Compar- 
ing the measured N(N n)/AT(N i), Ar(Alm)/Ar(Aln), and 
Af(Silll)/AT(Sill) ratios in Table 1] (all these ratios are < 1) 
with the values computed by CLOUDY for different densi- 
ties, we find that in this DLA logn(H)/cm~ 3 ^ —2 so 
that the ionisation corrections for carbon and nitrogen are 
|IC(C),IC(N)| < 0.1. This is also the case for the DLA 
in Q0913+072, although in this case we rely solely on the 
iV(Alm) /A r (Aln) ratio. Peroux et al. (2007) and Dessauges- 
Zavadsky et al. (2003) reached similar conclusions for the 
two sub-DLAs included in Table 1101 based on similar pho- 
toionisation analyses. For the two DLAs in the SDSS QSOs 
we have no comparable information and therefore cannot 
come to definite conclusions concerning ionisation correc- 
tions. However, even in these two cases, it is unlikely that 
to |IC(C), IC(N)| > 0.2. Larger ionisation corrections would 
imply volume densities n(H) < 0.001 cm -3 and therefore lin- 
ear sizes for the absorbing structures in excess of ~ 100 kpc 
(for JV(Hi) — 2 x 10 20 cm -2 ); it seems unlikely that gas in 
such large volumes would maintain internal velocity disper- 
sions of only a few km s _1 as required by the narrow line 
widths. 

In summary, our best estimates of the ionisation struc- 
ture of the sample of six DLAs and sub-DLAs consid- 
ered here leads us to conclude that the measured ratios 
AT(Cll)/iV(Ol) and iV(N i) /N(0 i) reflect the true values of 
C/O and N/O in the DLAs to within 0.2 dex and probably 
to within 0.1 dex in most cases. 



5 THE C/O RATIO AT LOW METALLICITIES 

One of the aims of the present work is to investigate the be- 
haviour of the C/O ratio at low metallicities, motivated by 
the results of Akerman et al. (2004) and Spite et al. (2005) 
which we now summarise. Akerman et al. (2004) measured 
the abundances of C and O in 34 halo stars of metallicities 
from [Fe/H] = —0.7 to —3.2 and combined such data with 
analogous measurements in disk stars from the literature. 
The pattern that emerged from their study is illustrated in 
Figure [7] Briefly, in halo stars [C/O] ~ —0.5 (correspond- 
ing to an observed log (C/O) ~ —0.77 since in the Sun 
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Figure 5. Ionisation corrections (defined as in eq. [TJ vs. gas 
density for observed ions which are dominant in H I regions. The 
values shown here are for the 2 abs = 2.07624 DLA in Q2206-199, 
but similar corrections apply to the other DLAs in Table ITOl 

log(C/O) = -0.27; Asplund et al. 2005) □ At the higher 
metallicities of most disk stars, C/O rises to solar propor- 
tions; in the Galactic chemical evolution models considered 
by Akerman et al. (2004) this rise is interpreted as the addi- 
tional contribution to carbon enrichment of the interstellar 
medium from the winds of Wolf-Rayet stars — whose mass 
loss rates are known to increase with metallicity — and, to a 
lesser extent, from the delayed evolution of low and inter- 
mediate mass stars. 

More puzzling is the apparent upturn in [C / O] with de- 
creasing metallicity when [O/H] ,$ — 1, since conventional 
stellar yields would predict the opposite behaviour: the ex- 
pected time lag in the production of C relative to O in the 
first episodes of star formation should cause the [C/O] ratio 
to decrease dramatically as we move to the lowest metal- 
licities in Figure [7] (interested readers are referred to the 
extensive discussion of the relevant chemical evolution mod- 
els by Akerman et al. 2004). Akerman et al. speculated that 
the trend towards higher, rather than lower, values of [C/O] 
when [O/H] < — 1 may be due to remaining traces of high 
carbon production by the first stars to form in the halo 
of the Milky Way proto-galaxy. Indeed, some calculations 
of nucleosynthesis by Population III stars (e.g. Chieffi & 
Limongi 2002) do entertain high carbon yields (relative to 
oxygen) that can reproduce the observations in Figure [7] 
More recently, Chiappini et al. (2006) have proposed that 
the observed behaviour can also be explained if stars of lower 
metallicity rotate faster: the higher rotation speeds greatly 



7 Note that in the earlier compilation by Grevesse Sz Sauval 
(1998) the solar ratio of C and O had a similar value log (C/0)g = 
—0.31; thus the C/O ratio appears to be robust the 3D and non- 
LTE effects which have introduced significant changes to the in- 
dividual abundances of C and O. 



Figure 6. Ratios of successive ion stages vs. gas density for the 
DLA in line to Q2206— 199. The values shown here were calcu- 
lated with a grid of CLOUDY models (see text for details). The 
measured values of the ion ratios (given in Table |9j indicate 
logra(H)/cm~ 3 ^ -2. 

increase the yields of C (and N) by the massive stars which 
are the main source of O. 

Some uncertainty in the interpretation of the results in 
Figure[7]is introduced by the fact that Akerman et al. (2004) 
based their abundance measurements on one-dimensional, 
local thermodynamical equilibrium (LTE), stellar atmo- 
sphere models. The authors pointed out themselves the pos- 
sibility that the neglect of non-LTE corrections to the high 
excitation C I and O I lines used for abundance determina- 
tions may mimic (or at least enhance) the rise in [C/O] with 
decreasing [O/H]. An initial investigation of such effects by 
Fabbian et al. (2006) seemed to dilute the significance of the 
results by Akerman et al. (2004) by showing that non-LTE 
corrections to the C I lines can be as high as —0.4 dex in stars 
of [Fe/H] ~ —3. A more recent study by the same group, 
however, has uncovered even larger corrections to the abun- 
dance of O at very low metallicities (Fabbian et al. 2008a), 
so that the case for a rise in [C / O] when [O /H] ^ — 1 is ac- 
tually reinforced when non-LTE effects on the formation of 
the stellar C I and O I lines are taken into account (Fabbian 
et al. 2008b). 

This discussion highlights the important role that low 
metallicity DLAs can play in clarifying the early nucleosyn- 
thetic sources of carbon. The physics of absorption line for- 
mation in diffuse interstellar clouds is considerably simpler 
than that pertaining to stellar atmospheres; we do not have 
to take into account the geometry of the absorbers nor their 
exact thermodynamical conditions. The only two limiting 
factors to the accuracy with which we can measure the C/O 
ratio are: (i) line saturation and (ii) ionisation corrections. 
In the preceding sections we have shown that the former, 
if present at all in the data considered here, would lead us 
to underestimate the value of C/O, while the latter would 
work in the opposite direction but the effect is unlikely to 
be greater than 0.1 dex for the DLAs in the present study. 
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Figure 7. [C/O] vs. [O/H] for Galactic stars from the study 
by Akerman et al. (2004). Black dots: disk stars; blue squares: 
halo stars. The green line shows the evolution of the [C/O] ra- 
tio with metallicity for a 'standard' Galactic chemical evolution 
model using conventional stellar yields (see Akerman et al. 2004 
for further details). A general feature of this class of models is 
the drop in [C/O] at the lowest metallicities, as the production 
of C in massive stars lags behind that of O. The stellar measure- 
ments by Akerman et al. (2004) provided the first indication of a 
different behaviour, suggesting that the [C/O] ratio may in fact 
be closer to solar at metallicities [O/H] ^ — 2.5. The disconti- 
nuity at [O/H] = —0.5 in the line describing the model evolution 
corresponds to the transition from the halo to the disk. 



When the values of [C/O] measured in the six 
DLAs/sub-DLAs in the present sample are compared with 
those of Galactic stars as in Figure [SJ it appears that both 
sets of data paint a consistent pictured The DLA measure- 
ments corroborate the conclusion by Akerman et al. (2004) 
in showing relatively high values of C/O in the low metal- 
licity regime and complement the still very limited statistics 
of this ratio in stars of metallicity [O/H] < —2. 

Although the high redshift DLAs and the halo stars ap- 
pear to constitute a continuous sequence in Figure [8] we 
warn against such a simplistic interpretation. The chemical 
enrichment histories of the Milky Way halo and of the galax- 
ies giving rise to DLAs need not be the same. Therefore it 
may well be that, as chemical evolution progresses, different 
DLAs follow different paths in the [C/O] vs. [O/H] plane, 
some similar to that of the Milky Way stellar populations 
while others may be different. Such a variety of behaviours 
has already been seen (in other elements) in Local Group 
galaxies (e.g. Pritzl et al. 2005). What our data show is that 
at metallicities below 1/100 of solar DLAs and Galactic stars 
concur in showing elevated ratios of C/O. Akerman et al. 
(2004) speculated that an approximately solar value of C/O 
might be recovered when [O/H] 5,-3; the new data pre- 
sented here are certainly consistent with this extrapolation 
of the Galactic measurements. 



8 The values of [C/O] plotted in Figure[8]were deduced from the 
values of log (C/O) listed in Table [TOl after subtracting the solar 
value log (C/0)q = —0.27 from the compilation by Asplund et al. 
(2005). The errors shown were obtained by combining in quadra- 
ture the la errors in the column densities iV(Cll) and W(Ol) 
returned by VPFIT and listed in Tables EJ [5j [T] and [9] or given 
by the original authors for the two DLAs from the literature. In 
all cases, a minimum column density uncertainty of 0.05 dex was 
assumed. 



Figure 8. The values of [C/O] deduced here for DLAs (red filled 
triangles) and sub-DLAs (open triangles) broadly agree with the 
stellar abundances measured at similar metallicities, supporting 
the suggestion by Akerman et al. (2004) that carbon and oxygen 
may be produced in near-solar proportions in the earliest stages of 
galactic chemical evolution. Other symbols have the same mean- 
ing as in Figure [7] 

6 THE N/O RATIO AT LOW METALLICITIES 

In contrast with the paucity of carbon measurements, sev- 
eral studies have targeted the abundance of N in DLAs since 
the first determination of the N/O ratio in a metal-poor 
DLA by Pettini, Lipman, & Hunstead (1995). These au- 
thors pointed out the potential of such data for clarifying 
the nucleosynthetic origin of N and unravelling the chemi- 
cal enrichment history of the galaxies giving rise to damped 
systems. Henry & Prochaska (2007) and Petitjean, Ledoux, 
& Srianand (2008) provide recent discussions of these issues 
and comprehensive reviews of earlier literature. 

At present, there are 35 DLAs with reliable measure- 
ments (or useful upper limits) of N/O, collected in Ta- 
ble lllrl Despite this body of data, it is still instructive to 
consider the three measurements of N/O reported here, as 
they refer to some of the most metal-poor DLAs known. 
In order to understand their importance, it is necessary to 
summarise current ideas on the production of N in stars of 
different masses — we refer the interested reader to Henry 
& Prochaska (2007), Pettini et al. (2002a), and references 
therein for more detailed treatments. 

The key points are as follows. There is general agree- 
ment that the main pathway for the production of N in 
stellar interiors is a six step process in the CN branch of the 
CNO cycle which takes place in the stellar H burning layer, 
with the net result that 14 N is synthesised from 12 C and le O. 
Less well established, however, are the N yields of stars of 
different masses and metallicities. The issue is complicated 
by the fact that N can be of either 'primary' or 'secondary' 
origin, depending on whether the seed carbon and oxygen 
are those manufactured by the star during helium burning, 
or were already present when the star first condensed out of 
the interstellar medium (ISM). 

This dual nature is revealed when the abundances of 

9 In most of these cases S is used as a proxy for O because, unless 
far-UV O I lines are accessible (and these are often blended in the 
Lya forest), the strong OlA1302 line is saturated in most DLAs. 
The use of S as a proxy for O is justified by the fact that S behaves 
like a typical alpha-capture element in Galactic stars (e.g. Nissen 
et al. 2007). 
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Table 11. N, O, S, and Fe Abundance Measurements in DLAs 



QSO 


2 abs 


log iV(Hi) 
(cm" 2 ) 


log JV(N i) 
(cm -2 ) 


logJV(Oi) 
(cm -2 ) 


logJV(Sn) 
(cm -2 ) 


log(0/H) + 12 a 


io g (N/or 


log (Fe/0) a 


Ref. b 


onono— 2fi2n 


3.3901 


21.41 


14.73 


16.45 




7.04 


—1.72 


— 1.58 


x 


om oo+i 300 


2 30903 


21.37 


15.03 




15.09 


7.22 


— 1.56 


—1.50 


2 


Om 1 2 — 306 


2.41844 


20.50 


13.16 


14.95 




6.45 


— 1.79 


— 1.62 


3 


0(1901 -U1 1 90 


3 38fi^Q 


21.26 


15.33 




15.21 


7 45 


— 1.38 


— 1.36 


4 


T0307— 4945 




20.67 


13.57 


15.91 




7.24 


—2.34 


—1.70 


5 


00347—381 Q 


3 02485 


20.73 


14.64 


16.45 


14.74 


7.72 


— 1.81 


—2.10 


3 


00407—441 fl 


2 5505 


21.13 


14.55 




14.82 


7.19 


— 1.77 


— 1.37 


g 


00407—441 D 


2.5950 


21.09 


15.07 




15.19 


7.60 


— 1.62 


— 1.54 


5 


00528—2505 


2.1410 


20.95 


14.58 




14.83 


7.38 


— 1.75 


— 1.48 


7 


HSD741 -U4741 


3.01740 


20.48 


13.98 




14.01 


7.03 


— 1.53 


— 1.46 


g 


0841 4-1 256 


2 37452 


20.99 


14.60 




14.69 


7.20 


— 1.59 


—1.43 


g 


f)£41 4-1 9^fi 


2 47621 


90 78 


1 3 Q4 


1 fi 1 ^ 


14 48 


7 ^7 
1 .O f 


—2 21 




n 
■ ' 




3.2458 


20.30 


14.15 




14.65 


7.85 


—2.00 


— 1.61 


g 




9 fS1 843 


20.34 


12.29 


14.63 




6.29 


—2.34 


— 1.64 


10 


Vc^ u o u 6 o m o 


3 2353 


20.18 


13.82 




13.67 


6.99 


— 1.35 




11 


Ol 1 PI8— D77 


3 60767 


20.37 


< 12.84 


15.37 




7.00 


< —2.53 


—1.49 


3 


Ol 21 (14-1 7 


1.89177 


20.63 


14.71 




14.96 


7.83 


— 1.75 


—1.45 


9 


01 9^94-nsi ^ 


9 33779 


20.80 


14.63 




14.83 


7.53 


— 1.70 


— 1.62 


7 


Ol H^l 4-1 7fl 


1.77637 


21.14 


< 15.23 




15.08 


7.44 


< — 1 35 


— 1.95 


2 


Ol 3374-1 1 3 


2.79581 


21 00 


13.99 


15 74 




6.74 


— 1.75 c 


— 1.41 


3 


Ol 4094-0Q c > 


9 4^690 


20.54 


< 13.19 


15.15 




6.61 


< —1.96 


—1.41 


12 


11 435+5359 


2.3427 


21.05 


14.67 




14.78 


7.23 


— 1.61 




s. 


J1443+2724 


4.224 


20.95 


15.52 




15.52 


8.07 


-1.50 


-1.69 


13 


SDSS J1558-0031 


2.70262 


20.67 


14.46 


15.87 


14.07 


7.20 


-1.41 


-1.76 


11 


SDSS J1558+4053 


2.55332 


20.30 


12.66 


14.54 




6.24 


-1.88 


-1.47 


10 


GB1759+7539 


2.62530 


20.76 


15.11 




15.24 


7.98 


-1.63 


-1.66 


8 


Q2059-360 


3.08293 


20.98 


13.95 


16.09 


14.41 


7.11 


-2.14 


-1.61 


3 


Q2206-199 


2.07624 


20.43 


12.79 


15.05 




6.62 


-2.26 


-1.72 


10 


Q2230+02 


1.8644 


20.85 


15.02 




15.29 


7.94 


-1.77 


-1.30 


9, 8 


Q2231-002 


2.06616 


20.53 


< 15.02 




15.10 


8.07 


< -1.58 


-1.77 


2 


HE2243-6031 


2.33000 


20.67 


14.88 




15.02 


7.85 


-1.64 


-1.60 


15 


Q2332-094 


3.05723 


20.50 


13.73 


15.95 


14.37 


7.45 


-1.94 c 


-1.58 


3 


Q2342+342 


2.9082 


21.10 


14.92 




15.19 


7.59 


-1.77 


-1.67 


8 


Q2343+1232 


2.43125 


20.35 


15.22 




14.86 


8.01 


-1.14 


-1.70 


16 


Q2348-1444 


2.27939 


20.59 


13.35 




13.75 


6.66 


-1.90 


-1.41 


9 



a When the oxygen abundance is not available, S has been used as a proxy for O by assuming the solar ratio log(0/S)g = +1.50 
(Asplund et al. 2005). 

b References — 1: Molaro et al. (2000); 2: Dessauges-Zavadsky et al. (2004); 3: Petitjean, Ledoux, & Srianand (2008); 4: Ellison ct al. 
(2001); 5: Dessauges-Zavadsky et al. (2001); 6: Lopez & Ellison (2003); 7. Centurion et al. (2003); 8. Henry & Prochaska (2007). 9. 
Dessauges-Zavadsky et al. (2006); 10: This work; 11. Lu, Sargent, & Barlow (1998); 12. Pettini et al. (2002a). 13. Ledoux et al. (2006b). 
14. O'Meara et al. (2006). 15: Lopez et al. (2002); 16: D'Odorico, Petitjean, & Cristiani (2002). 
c (N/0) ratio only for components in common between Ni and Ol. 



N and O in extragalactic Hn regions (these are the main 
source of measurements of these two elements in the nearby 
universe) are compared as in Figure [9F°I When the oxy- 
gen abundance is greater than about half solar, that is 
when log(0/H) + 12 > 8.3 adopting log(0/H) Q + 12 = 
8.66 from Asplund et al. (2005), N/O rises steeply with 
increasing O/H; this is the regime where N is predomi- 



10 See Pettini et al. (2002a) for the sources of the Hn region 
measurements reproduced in Figure [9] Although much more ex- 
tensive data sets have become available recently through large 
scale surveys such as the Sloan Digital Sky Survey (e.g. Mallery 
et al. 2007; Liang et al. 2007 — see also Nava et al. 2006), the locus 
occupied by Hn regions in the log (N/O) vs. log (O/H) plane has 
not changed significantly. 



nantly secondary!"] At lower metallicities on the other hand 
(log (O/H) + 12 ^ 8.0), N is mostly primary and its abun- 
dance tracks that of O; this results in the flat plateau at 
log (N/O) ~ -1.5 evident in Figure [9] 

Primary N is thought to be synthesised most effectively 
by intermediate mass stars on the asymptotic giant branch; 
consequently, its release into the interstellar medium pre- 
sumably takes place some time after the massive stars which 
are the main source of oxygen have exploded as Type II su- 
pernovae. Henry et al. (2000) calculated the delay to be 

11 Henry, Edmunds, & Koppen (2000) pointed out that the rise 
in N/O with O /H is steeper than would be normally expected for 
a purely secondary element and proposed that it is augmented by 
a decreasing O yield with increasing metallicity. 
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Figure 9. Abundances of N and O in extragalactic Hn regions (small dots) and high redshift damped Lya systems (triangles). The 
former are from the sources assembled by Pettini et al. (2002a), with the addition of more recent data from Melbourne et al. (2004). 
The DLA values are from the compilation in Table 1111 filled triangles denote DLAs where the abundance of O could be measured 
directly, while open triangles are cases where S was used as a proxy for O. The three green filled triangles are the new determinations 
reported in this paper. The error bars in the bottom right-hand corner give an indication of the typical uncertainties; the large yellow 
dot corresponds to the solar abundances of N and O (Asplund et al. 2005). The dashed lines are approximate representations of the 
primary and secondary levels of N production (see text). 



At ~ 250 Myr for a 'standard' stellar initial mass function 
(IMF) and published stellar yields. 

Damped Lya systems add information that can be of 
value in assessing the validity of this overall picture. First 
of all, they offer measures of the abundances of O and N 
in entirely different environments from nearby star-forming 
galaxies. Second, they sample low metallicity regimes which 
are rare today. And third, by virtue of their being at high 
redshift, they can in principle provide a better sampling of 
the time delay in the production of primary N relative to 
O — the canonical At ~ 250 Myr is a much larger fraction of 
the time available for star formation at z = 2 — 3 than today. 
In other words, the chances of 'catching' a galaxy in the 
interim period following a burst of star formation, when the 
O from Type II supernovae has already dispersed into the 
interstellar medium but lower mass stars have yet to release 
their primary N, are considerably higher in DLAs than in 
metal-poor star-forming galaxies in the local universe. Such 
a situation should manifest itself as a displacement below 
the primary plateau in the log (N/O) vs. log (O/H) plane 
in Figure 

The number of DLAs with reliable measures of the N 
and O abundances has slowly grown since the first deter- 
mination by Pettini et al. (1995); the current total of 35 
such systems (Table lll[) provides a moderate size sample 
with which to confront the above scenario. When DLAs are 
compared with present-day galaxies, as in Figure [5] it can 
be seen that in most cases their values of N/O fall close to 
the plateau defined by metal-poor H II regions. If primary N 
and O are synthesised by stars of different masses, as gener- 
ally believed, the good agreement between local and distant 
galaxies attests to the universality of the IMF — a top-heavy 



IMF in DLAs would result in lower values of N/O, contrary 
to observations. 

Out of the 35 DLAs in our sample, 12 (~ 34%) lie more 
than 0.3 dex [twice the typical la error in the determina- 
tion of log (N/O)] below the primary plateau. For compar- 
ison, At/[1/H(z)\ ~ 0.18 at the mean (z) = 2.70 of the 35 
DLAs in Table [TT1 where 1/H(z) is the age of the universe 
at redshift z (in today's consensus cosmology). The differ- 
ence between these two fractions is probably not significant, 
given that the DLA sample in Table [TT1 is somewhat biased 
towards chemically unevolved DLAs and that the age of the 
universe is in any case an upper limit to the time available 
for star formation. Certainly, the current DLA statistics are 
consistent with a time delay At of a few 10 s years and allay 
earlier concerns based on more limited samples (e.g. Pettini 
et al. 2002a). 

Figure [9] also emphasises the importance of extending 
measurements of N/O to the lowest metallicity DLAs: as the 
gap between secondary and primary N production widens, 
it becomes easier to recognise additional sources of primary 
nitrogen, if they exist. That is, if intermediate mass stars 
are the sole producers of primary N, it should be possible to 
find DLAs with values of N/O close to the secondary line in 
Figure On the other hand, if massive stars were able to 
synthesise some primary N, we might expect to see a mini- 
mum value of N/O, corresponding to the N ejected into the 
ISM at the same time as O. This idea was first proposed 
by Centurion et al. (2003 — see also Molaro 2003) and re- 
cent work on stellar yields including the effects of rotation 
(e.g. Hirschi et al. 2007) suggests that it may be plausible. 
In this respect, the new data presented here (shown by the 
green triangles in Figure [9} add some weight to the sug- 
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gestion by Centurion et al. (2003): the lowest metallicity 
DLAs show values of N/O as low as, but no lower than, 
the lowest previously measured (in more metal-rich DLAs). 
While the statistics are still very limited, the available data 
are at least consistent with a 'floor' in the value of N/O at 
log (N/O) ~ —2.3. If this does indeed reflect the primary N 
yield by massive stars, it would imply that at these metal- 
licities massive stars can account for ~ 15% of the total 
primary N production. Clearly, it will be of interest to test 
the validity of this interpretation as more metal-poor DLAs 
are studied at high spectral resolution and the sample of 
N/O measurements at [O/H] < —2 increases in the years 
ahead. 



7 SUMMARY AND CONCLUSIONS 

In this paper we have presented new observations of four 
damped absorption systems (three DLAs and one sub-DLA 
according to the conventional division between the two at 
log AT (Hi) = 20.3) which are among the least metal-enriched 
known. We have focussed in particular on the clues which 
the abundance patterns of these DLAs provide on the nu- 
cleosynthesis of the C, N, O group at low metallicities. Our 
main conclusions are as follows. 

1. DLAs at the low end of the metallicity distribution can 
be preselected from medium resolution QSO spectra by tar- 
geting damped Lya lines which apparently have no associ- 
ated metal absorption lines. As expected, such systems turn 
out to have narrow metal lines when observed at higher spec- 
tral resolution, consisting of only a few components spread 
out over a few tens of km s _1 . 

2. The most metal-poor DLAs known have oxygen abun- 
dances of ~ 1/250 solar ([O/H] ~ —2.4), and 'metallicities' 
which are lower by factors of ~ 2 — 3 when measured in 
carbon and iron. This level of chemical enrichment is sim- 
ilar to that typical of intergalactic Lya forest clouds and 
suggest that these DLAs have experienced very little, if any, 
local star formation activity. No DLAs have yet been discov- 
ered with [Fe/H] < —3, comparable to the iron abundance 
of the most metal- deficient stars in the Galactic halo, but 
such stars are very rare and larger samples of DLAs may be 
required to intersect such pristine clouds of neutral gas at 
high redshift. 

3. Together with two sub-DLAs from earlier published 
studies, we have assembled a small sample of six absorption 
systems where the C/O ratio can be measured with confi- 
dence. We find good agreement between the relatively high 
values of C/O in these metal-poor DLAs and those mea- 
sured by Akerman et al. (2004) in halo stars of similar oxy- 
gen abundance. The DLA determinations are less ambigu- 
ous than their stellar counterparts. The agreement we have 
found suggests that corrections to the stellar abundances of 
C and O resulting from more sophisticated treatments of the 
stellar atmospheres (including 3-D and non-LTE effects) are 
likely to be of similar magnitude for the lines analysed by 
Akerman et al. (2004), so that the net corrections to the 
values of C/O they derived may not be large. The apparent 
rise of the C/O ratio as we move from stars with [O/H] ~ — 1 
to lower metallicities points to an additional source of car- 
bon from the massive stars responsible for early episodes of 



nucleosynthesis, in the Milky Way halo and in high redshift 
DLAs. 

4. The data presented here extend the growing sample of 
N/O measurements in DLAs to lower oxygen abundances, 
where the distinction between secondary and primary chan- 
nels widens. Our results strengthen earlier suggestions of a 
'floor' in the value of N/O at log (N/O) ~ -2.3 (or ~ 1/25 
of the solar N/O ratio), although larger samples of metal- 
poor DLAs are necessary to establish its reality. If this floor 
reflects the small primary production of N by massive stars, 
it amounts to ~ 15% of the total primary N yield integrated 
along the IMF. If most of the primary N is from intermediate 
mass stars, the generally good agreement in N/O between 
most DLAs and nearby low metallicity galaxies is an indica- 
tion that the IMF does not vary significantly between these 
different environments and cosmic epochs. More generally, 
our results add weight to the universality of the stellar yields 
of C, N, and O (Henry 2004), given the consistency in the 
ratios C/O and N/O between high redshift DLAs, Galactic 
stars, and nearby Hn regions. 

What is still unclear is where and when these trace met- 
als, which account for less than ~ 5 x 10~ 5 of the baryons 
by mass, were synthesised. The element ratios in the most 
metal-poor DLAs offer clues which may also be relevant to 
the origin of the metals polluting the intergalactic medium 
(IGM), where abundances are less well determined (as they 
rely on photoionisation models) but appear to be similar to 
those deduced here. 

In particular, it is now well established that star- 
forming galaxies at z — 2 — 3 drive large-scale gas outflows 
into their surroundings (Adelbeger et al. 2005; Steidel et al. 
in preparation) ; the metallicity of the outflowing material is 
likely to be within a factor of a few of solar (Pettini et al. 
2002b; Erb et al. 2006a). A novel suggestion recently put 
forward by Schaye, Carswell, & Kim (2007) is that the low 
metallicities commonly deduced for the IGM are mislead- 
ing, resulting from the poor mixing between small parcels of 
ejected gas of near-solar composition and true intergalactic 
gas of primordial origin on much larger scales. The relatively 
high C/O values we have uncovered in the DLAs studied 
here could be viewed in this context, as they are not dissim- 
ilar to those also measured in disk stars with [O/H] £ — 0.5 
(see Figure [8} . In this case, the good agreement with metal- 
poor halo stars with [O/H] £ — 2 would be coincidental. A 
more extensive set of C/O measurements in both stars and 
DLAs should clarify the situation; for example, a definite 
trend of increasing C/O with decreasing [O/H] would not 
be expected on the basis of the model proposed by Schaye 
et al. (2007). 

The few DLAs with the lowest values of N/O are more 
discriminatory in this respect. The median ages of the galax- 
ies which support the outflows at z — 2 — 3 are between 
~ 300 and ~ 600 Myr (Shapley et al. 2001; Papovich et 
al. 2001; Erb et al. 2006b), comparable to, or somewhat 
larger than, the time delay At for the release of primary 
N from intermediate mass stars. But even in the younger 
galaxies, the secondary N produced when the metallicity is 
high ([O/H] > -0.5, Erb et al. 2006a) would result in N/O 
ratios significantly higher than observed in the most metal- 
poor DLAs. For example, in the young Lyman break galaxy 
MS1512-cB58, where [0/H]~ -0.4 (Teplitz et al. 2000), 
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Pettini et al. (2002b) deduced log (N/O) = -1.8, a factor 
of ~ 3 higher than the 'floor' at log (N/O) ~ —2.3 exhibited 
by DLAs. Taken together, the C/O and N/O measurements 
reported here favour an origin of these elements in stars of 
genuinely low metallicity, rather than in diluted ejecta from 
galaxies which have already attained high levels of chemical 
enrichment. 

In this case, our data hint to the possibility that massive 
stars may contribute significantly to the nucleosynthesis of 
both C and N, as well as O, at the lowest metallicities. Stellar 
rotation may be the common factor here, through the mixing 
it induces in the stars' interiors, but it remains to be estab- 
lished if physically motivated metallicity-dependent C, N, 
and O yields can reproduce the relative abundances of these 
elements measured in DLAs, when incorporated into ap- 
propriate chemical evolution models. From an observational 
point of view, it should be possible with current instrumen- 
tation and some dedicated effort to increase the numbers of 
absorption systems with reliable measures of C/O and N/O 
and place the conclusions reached here on stronger statis- 
tical footing. We look forward to such developments in the 
years ahead. 
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